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We made deep near-infrared (JHKs) imaging polarimetry toward the Ser- 
pens cloud core, which is a nearby, active cluster forming region. The polar- 
ization vector maps show that the near-infrared reflection light in this region 
'p_i! mainly originates from SVS2 and SVS20, and enable us to detect 24 small in- 

frared reflection nebulae associated with YSOs. Polarization measurements of 
near-infrared point sources indicate an hourglass-shaped magnetic field, of which 
I symmetry axis is nearly perpendicular to the elongation of the C^®0 ( J = 1 — 0) 

or submillimeter continuum emission. The bright part of C^^O (J = 1 — 0), 
submillimeter continuum cores as well as many class O/I objects are located just 



O 



> 

G> I toward the constriction region of the hourglass-shaped magnetic field. Applying 

^ I the Chandrasekhar & Fermi method and taking into account the recent study on 

I the signal integration effect for the dispersion component of the magnetic field, 

^ I the magnetic field strength was estimated to be ~100 /iG, suggesting that the 

O I ambient region of the Serpens cloud core is moderately magnetically supercritical. 
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These suggest that the Serpens cloud core first contracted along the magnetic 
field to be an elongated cloud, which is perpendicular to the magnetic field, and 
that then the central part contracted cross the magnetic field due to the high 
density in the central region of the cloud core, where star formation is actively 
continuing. 

Comparison of this magnetic field with the previous observations of molecular 
gas and large-scale outfiows suggests a possibihty that the cloud dynamics is 
controlled by the magnetic field, protostellar outfiows and gravitational infiows. 
Furthermore, the outfiow energy injection rate appears to be larger than the 
dissipation rate of the turbulent energy in this cloud, indicating that the outfiows 
are the main source of turbulence and that the magnetic field plays an important 
role both in allowing the outfiow energy to escape from the central region of the 
cloud core and enabling the gravitational infiows from the ambient region to the 
central region. These characteristics appear to be in good agreement with the 
outfiow-driven turbulence model and imply the importance of the magnetic field 
to continuous star formation in the center region of the cluster forming region. 

Subject headings: circumstellar matter infrared: stars ISM: individual (Serpens) 
ISM: magnetic fields polarization stars: formation 



1. Introduction 



Stars are formed by gravitation in molecular clouds having both turbulence and magnetic 
fields in the Galaxy, an d most of stars are thought to be formed in clusters (e.g.. lLada fc Lada 
20031 : lAllen et al.ll2007l ). A mass spectrum of prestellar condensations is repo rted to have the 
powe r similar to that of the stellar IMF both in dust continuum obser yations (IReid fc Wilson 



20061 . and references therein) and molecular-line observations (e.g., Ilkeda et al.l 120071 ). and 



theoretical studies of turbulent molecular clouds ( iKlessen et al.lll998l and subsequent works) 
suggest that these condensations were formed through turbulent shock. One of the most 
promising sources of ordinary turbulence is outfiows from protostars, which are ubiquitous 
in star forming regions and are believed to be formed through the mediation of magnetic 
field. Magnetic fields are also considered to play an important role in dynamical evolution 
of molecular clouds and control of star formation, i.e., formation of molecular cloud cores 



and their collapse (e.g., iMcKee fc Ostrikerl 120071 ). 



Recently, iLi &: Nakamural (|2006[ ) and iNakamura &: Li I (120071 ) presented realistic 3D 
MHD simulations of cluster formation, taking into account the effect of protostellar outflows 
as well as initial turbulence and a magnetic fleld. In their simulations, they indicated that 
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the initial turbulence is quickly replaced by turbulence generated by protostellar outflows, 
keeping the quasi-equilibrium state with a slow rate of star formation, and that magnetic 
fields are dynamically important if their initial strengths are not far below the critical value 
for static cloud support because of the amplification by the outfiow-driven turbulent motions. 
The magnetic field is expected to infiuence the directions of outfiow ejection and propagation 
and the transmission of outflow energy and momentum to the ambient medium. However, 
the magnetic fleld structures have not always been observationally clear in/around cluster 
forming regions, particularly around nearby cluster forming regions because of the lack of 
deep, wide-fleld near-infrared (NIR) polarimetry data. 

The Serpens cloud core is one of the nearbj0, active low-mass star forming regions 
at the northern p art of the Serpens cloud and many observational work s have been done 



(lEiroa et al. 



Harvey et al. 



20081. and references there in). Recent mid-IR studies (e.g.. iKaas et al.l 12004 



20071 : IWinston et al.ll2007l ) revealed that a lot of embedded young stellar objects 



(YSOs), including Class O/ I objects, are located toward an aggregate of (sub) m illimeter 



dust continuum cores (e.g., iDavis et al.l Il999l : iKaas et al.l 120041: Enoch et al.l 120071 ) . which 



Olmi fc Testi 



2002 



consists of two sub-clumps (NW and SE sub-clumps; 
region and is enveloped by ambient molecular gas (e.g. ^^CO, and C^^O 
2000l : loimi fc Testill2002l ). 



in the central 



McMulhn et al. 



Many outflow activities that are related to star formation have been taking place in the 
Serpens cl oud core. CO high velocity flows are reported to be widely spread over the cloud 



core (e.g., iWhite et al.lll995l : iDavis et al.lll999l : iNarayanan et al.ll2002l ). Compact molecular 
outflows of higher density tracers and H2 jet-like knots are associated with the submil- 



limeter cores (e.g.. 


Curiel et al. 


1999; 


Williams & Mvers 


2000 





Curiel et al.l Il996l : iHerbst et al.l Il997l : IWolf-Chae et al.l Il998l : iHodapp 



PA~ 155° on an average with deviation of a few 10° (see Table 5 of lOlmi fc Testil |2002| ). 



whic h is nearly parallel to the alignment direction, from NW to SE, of the two sub-clumps 



Davis et al.l Il999t iKaas et al.l l2004j : Enoch et al.l l2007l). or to the cloud elongation in 



le.g., 

^^CO, C^^O, and o ther higher d ensity tracers 



e^ 



McMuUin et al. 2000; Olmi fc Testi 



2OO2I ). In addition, IPavis et al.l ^9M) and Iziener fc Eisloffell jl999 [rfound. through opti 



cal narrow-band imaging, that many HH objects emanate from the two sub-clumps to the 
ambient region, penetrating the dense part of the central region. 



The Serpens Reflection Nebula (SRN) illuminated by SVS 2 ( IStrom et al.lll976l ) has been 



^We assume a distance of ^260 pc for the Serpens cloud , follow ing the most of the recent papers on the 
Serpens cloud and based on the discussion of IStraizys et al.l (|2003l ) on the center distance of the Aquila Rift 
system. 
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extensively s tudied by polarimetric measurements both in optical and near-infrared 



1997; Huard et al 



wave l engths (iKing et aUl983l : iWarren-Smith et al.lll987l : iGomez de Castro et al 



19971 ). NIR polarimetric measurements ( ISogawa et al 



1997 



1988 



(NIR) 



Sogawa et al 



Huard et al. 



Gomez de Castro et al 



19971) probed also some other obscured reflection nebulae around SVS 2 in detail 
( jl988 ) suggested the magnetic field of a NW-SE direction based on the elongation of the 
refiection nebulae around several YSOs in the central region of Serpens cloud core. In con- 
trast, the NIR polarization measurement of a background star candidate suggested rather 
different direction of ma gnetic field because its polarization angle was nearly perpendicular 
to the NW-SE direction (jSogawa et al.l 119971 ). However, this measurement was only for one 
background candidate, which in fact has a possibility of being a YSO in the Serpens cloud 
core and its polarization originating from the YSO itself. Therefore, it is vitally important 
to measure more background stars to resolve this discrepancy and to know the magnetic 
field structure toward the Serpens cloud core. 

We conducted deep, JHKs imaging polarimetry of the Serpens cloud core to reveal the 
magnetic field structure in this region. We also searched for more NIR refiection nebulae 
associated with YSOs. Here, we present the results of our imaging polarimetry in the Serpens 
cloud core by comparing the data from the previous observations and discuss the role of the 
magnetic field in this region. 



2. Observations and Data Reductions 

Toward the Serpens cloud core (Figured]), simultaneous JHKs polarimetr ic observations 



were carried out on 2006 June 23 UT with the imaging polarimeter SIRPOL (IKandori et al. 



20061 ). which is an attachment of the near-infrared camera SIRIUS mounted on the IRSF 1.4- 
m telescope at the South Africa Astronomical Observatory. The SIRIUS camera is equipped 
with three 1024 x 1024 HgCdTe (HAWAII) arr ays, JHKs filters , and dichroic mirrors. 



which enables simultaneous JHKs observations (iNagashima et al.l Il999l : iNagayama et al. 



20031 ) . The field of view at each band is ~7f7 x 7'7 with a pixel scale of 0'.'45 pixel 
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We obtained 10 dithered exposures, each 10 s long, at four wave-plate angles (0°, 22.5°, 
45°, and 67.5° in the instrumental coordinate system) as one set of observations and repeated 
this 9 times. Sky images were also obtained in between target observations. Thus, the total 
on-target exposure time was 900 s per wave-plate angle. The seeing was ~ 1.2" at i^s during 
the observations. Twilight fiat-field images were obtained at the beginning and end of the 
observations. 



Standard procedures, dark subtraction, fiat-fielding with twilight-fiats, bad-pixel sub- 
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stitution, sky subtraction, and averaging of dithered images, were applied with IRAF. 
We first calculated the Stokes parameters as follows; Q = Jq — 145, U = I22.5 — hi. 5, 
/ = (/o + I22.5 + hb + -^67.5)/2, where Jq, 122.55 hb, and /67.5 are intensities at four wave-plate 
angle s. To obtain the Sto kes parameters in the equatorial coordinate system, a rotation of 
105° (IKandori et al.ll2006l ) was applied to them. We calculated the degree of polarization 
P, and the polarization angle 9 as follows; P = ^/Q^-TLP/I, 6 = (1/2) tan-\U/Q). The 
polarization intensity (PI) is obtained by multiplying the total intensity (/) by the degree 
of polarization (P). The absolute accuracy of the position angle of polarization was est i- 
mated to be better than 3° at the first light observation of SIRPOL (IKandori et al.ll2006l ). 
The polarization efficiencies are 95.5%, 96.3%, and 98.5% at J, H, and Ks, respectively, 
and the instrumenta l polarization is less than 0.3% all over the field of view at each band 
( IKandori et al.l 120061 ). Due to these high polarization efficiencies and low instrument polar- 
ization, no particular corrections were made here. 

Aperture polarimetry was performed for H and Ks band point sources detected by 
DAOFIND in the field of view. No polarimetry for J band sources was done due to their 
much smaller number, compared with those of H and Ks band sources (see Figure [2] a, c, and 
e). APHOT of the DAOPHOT package was used to evaluate the point source magnitudes 
for four wave-plate angles at H and Ks. An aperture radius of 3 pixels was adopted for 
each band. The errors of the degree of polarization {AP) and the position angle were 
calculated from the p hotometric errors, and th e degrees of polarization were debiased as 
Pdebias = VP^ - AP^ (IWardle fc Kronberglll974h . Hereafter, we use P as substitute for this 
debiased value for the aperture photometry data. Only the sources with ph otometric errors 



of <p .l mag and P/ AP > 2 were used for analysis. The 2MASS catalog (jSkrutskie et al. 



20061 ) was used for absolute photometric calibration. The limiting magnitudes at 0.1 mag 
error level were estimated to be 18.6 at H and 17.5 at Ks. 



3. Results 



3.1. Polarizations of extended emission 



The JHKs polarization vector maps of the Serpens cloud core are shown, superposed 
on the total and polarization intensity images in Figure [2l These polarization maps clearly 
indicate that the central part of the refiection nebula is illuminated mainly by two sources; 
the north part (SRN) by SVS 2 a nd the south part by SVS 20, at ij" and Ks with two 
centrosymmetric patterns (see also ISogawa et al.l Il997l : iHuard et al.l 1 19971). while at J only 
SRN is dominant (see also ISogawa et al.lll997l) as is seen in the optical (IWarren-Smith et al. 



1987; Gomez de Castro et al. 



19881 ). This invisibleness at shorter wavelengths suggests that 
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the southern part of SRN around SVS 20 is more obscured than the northern part a round 
SVS 2, consistent with the Ay map deduced from H — K color (IHuard et al.l 119971 ) . The 



centrosymmetric patterns are more clearly shown in Figure [31 which is a Ks polarization 
vector map shown with a resolution higher than Figure [2fe. 

The PI images and vector patterns of SVS 2 clearly show that SVS 2 is associated with 
a bipolar structure with a dark lane. In the JHKs intensity images of Figures |2] a, c, and 
e, at shorter wavelengths, the NW lobe of the bipolar nebula is brighter than the SE lobe, 
while at longer wavelengths the SE lobe is brighter. This suggests that the NW lobe is near 
side and that the SE lobe is far side. The nebula structure and dark lane of SRN have been 



already r eported in the two polarimet ri c stud ies ( ISogawa et al.lll997l : iHuard et al.lll997l ). In 
addition, iPontoppidan &: DuUemond I (|2005[ ) modeled SRN as a disk shadow system with 
their imaging data, and suggested that SVS 2 is associated with a small disk, which is not 
unresolved, and a spherically symmetric envelope. 

The nebula illuminated by SVS 20 is clearly recognized at H and Ks with a centrosym- 
metric pattern around SVS 20. This object has a peculiar morphology with a ring and two 
arms protruded from that ring. Because we plan to report its details in a separate paper, 
including other YSOs with reflection nebulae, we will not mention the details here. 



3.2. Polarizations of nebulosities associated with YSOs 

3.2.1. The central region 

Figure [3]presents an higher resolution Ks polarization vector map with 3x3 pixel binning 
toward the central region of the overall image, superposed on the Ks intensity map. With 
this map and/or the highest resolution vector maps without binning, we identified stellar 
sources having reflection nebulae locally illuminated by themselves with centrosymmetric 
or centrosymmetric-like patterns. It is not easy to identify such sources only from Figure 
[3] due to the strong contamination from SVS 2 and SVS 20. It is also not easy toward 
the SVS 4 cluster, which is a compact cluster located to the south of SVS 20, due to the 
source congestion. We therefore used the highest resolution vector maps without binning for 
the sources having the strong contamination (see Appendix). The identified sources with 
reflection nebulae illuminated by themselves are marked in Figure [3], including SVS 2 and 
SVS 20, and are listed in Tabled 

Most of the identified sources are relatively bright in the central region. This is probably 
due to the strong light from SVS 2 and SVS 20 and only brighter sources with reflection 
nebulae may be detected. Except EC117 (SSTc2dJ 18300065+0113402), all the identified 



- 7- 



sources are classified as sources that have outer disks with an excess at least 8 and /or 2 4 
/im, i.e., class O/I, flat spectrum, class II and transition disk sources (jWinston et al.l 120071 ). 
Although ECU 7 is classified as a class III source wi thout an outer disk due to no detection 
of 24/im continuum (Table 4 of IWinston et al.ll2007l). it was reported that EC117 has a flux 
of 3.20±1.43 mJy at 24/im (Table 3 of iHarvey et al.l 120071 ). This could suggest the outer 
disk of EC117, but the signal to noise ratio of ~2 is not high enough for the robust detection 
at 24/im. Almost all members of the SVS 4 cluster seem to be associated with reflection 
nebulae. 



3.2.2. The NW region 

Figure m presents a high resolution i^s polarization vector map without binning toward 
the NW region of SRN. Here we identified stellar sources associated with self-luminous neb- 
ulae as well as those with reflection nebulae by using this map and listed them in Table 
[2J Three sources, DEOS, EC53, and EC67, are associated with reflection nebulae having 
centrosjTiimetric or partially centrosymmetric vector patterns. The other sources are asso- 
ciated with elongated nebulae or jet-like knots emanating from the sources in a straight line 
and their polarization vectors are almost perpendicular to their elongation directions. The 
elongated structures or knots are likely to be created/excited by outflows from these sources. 

The jet-like knots are clearly seen near the north-west of EC41, which was considered 



to be an embedded star but not a driving source of this jet (lEiroa fc Casalilll989l : iHodapp 



9991). These i et-like knots are reported to be mostly H2 emission with weak continuum 
(jHodapp I I1999I ). and their polarization of these knots is nearly perpendicular to the jet 
elongation, though the polarization directions are more scattered in the northern knots than 
in the southern knots of this jet-like structure. At if-band, the polarization vectors of 
weak continuum emission of the southern knots are also nearly perpendicular to the jet 
elongation, which is parallel to the radial direction from SMMl-FIRSl, not from EC41. 
Thus, SMMl-FIRSl is the illuminating source of these jet-like knots, and the jet-like knots 
could correspond to the cavity walls that were created by the outflow from SMMl-FIRSl. 

The jet-like structure from SMMl-FIRSl seems to continue farther away to a bow-shock- 
like nebulosity located at ~80-90" north-west of SMMl-FIRSl (or at ~ 20" north-west of 
EC28). The polarization vectors at this bow nebulosity indicate either SMMl-FIRSl or EC28 
is illuminating or exciting it. No information is available on whether the bow neb ulosity is 



really shock-excited H2 emission or not, due to its position outside Figure 3 of iHodapp 



(I1999I ). although week ff-band continuum emission is detectable with polarization angles 



similar to those of /Ts-band. The alignment with the jet structure, the bow nebulosity and 
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HH 460, which is located at ~ 4' north-west to SMMl-FIRSl flDavis et al.lllQQQf ). gives a 
hint that the bow nebulosity is related to the outflow frorn SMMl-FIRSl. The associations 
of the blueshifted CO l obe with HH 460 (jPavis et al.lll999l ) and of the bow nebulosity with 
the CS emission (CSl; iTesti et al.ll2000l ). which is considered to be related to the outflow, 
support the shock excitation of the bow nebulosity. However, it is impossible to completely 
exclude the possibility that EC28, which is the closest NIR source to the bow nebulosity, 
or SMMl-FIRSl itself contributes to the illumination of the bow nebulosity, due to the 
scattering of the polarization vector directions. In the midway from SMMl-FIRSl to this 
bow nebulosity , there exist some faint knots that are almost H2 emission (see Figure 3 of 
Hodapp Ill999l ). but no polarization is detectable for these knots. 



Th e chain of nebulosity knots, located just sout h-east to the 3 mm continuum core 
S68Nc JTesti fc Sargentlll998l : Iwilhams fc Mverslboooh. was report ed to be H2 emission jets 
that originate from the 3mm subknot a3/S68Nc (IHodapp I Il999l ). Although only several 
polarization vectors are shown toward these knots, they could imply that their polarization 
direction is nearly perpendicular to the jet elongation. 

A nebulosity protruding from EC38/S68Nb is seen, and its polarization vectors appear 
to be nearly perpendicular to the protruding direction. A faint, small, elongated nebulosity is 
recognized just south-east to SMMIO-IR. Although some nebulosities illuminated from SVS 
2 are also seen near SMMIO-IR, this nebulosity is most likely a nebulosity related to SMMIO- 
IR based on its morphology. No information is available on whether these two nebulositie s 
are shock-excited H2 emission or not, because they are out of Figure 3 of iHodapp I (Il999l ). 
We note that no H-hand emission is detectable toward SMMIO-IR, while very week if-band 
continuum is seen toward the nebulosity protruding from EC38/S68Nb. 



3.3. Polarizations of point sources 

We have measured H and Ks polarization for point sources, in order to examine the mag- 
netic field structures. Only the sources with photometric errors of <0.1 mag and P/AP > 3 
were used for analysis. 

The top panel of Figure [5] presents the polarization degrees at H versus H — Ks color 



diagra m for sources having polarization errors of < 0.3%. YSOs identified by lWinston et al. 



(120071 ) are not included in this diagram. In this diagram, the maximum of polarization degree 
ai an H — K color is roughly proportional to the H — K color, i.e., the extinction, having 
consistency that the origin of the polarization is dichroic absorption. Therefore, we consider 
the polarization of these point sources as the polarization of the dichroic origin, and that 
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their polarization vectors represent the directions of the local magnetic field averaged over 
their line of sight of the sources. In the nearby star forming regions such as the Taurus and 
Ophiuchus clouds, the highest value of the maxir num polarizatioii effici ency was reported to 
be P(H)/E(H - K) = 4.6 or P ( H)/A (H) = 1.6 JKusakabe et al.lboosh . which were derived 
from the data of IWhittet et al.l (120081 ). In Figure [5], a dashed line represents P{H)/[H — 
Ks) = 4.6 where the offset of the intrinsic H — Ks color is ignored. Our sources have the 
maximum polarization efficiency of P{H)/{H — Ks) = 6.2 (thick line) similar to that of the 
nearby star forming regions, and this is also consistent with the dichroic origin. 

The bottom panel of Figure [5] shows the H-hand polarization angles of the point sources 
with P < 6.2{H — Ks), of which the polarization vectors are shown in Figure [61 YSOs are 
not included in the bottom panel, but included in Figure O 

The polarization angles are mostly in a range of ~0-140° and their median and average 
angles are 63.5° and 64.6±35.6°, respectively. While the polarization angles are largely 
scattered, there is a tendency that the degree of scatter becomes smaller in the redder 
H — Ks color region. This tendency suggests that the polarization angles are more confined 
in the inner region (redder color region) than in the outer region. 

The magnetic field is neither simply straight nor random over the whole field (Figure 
|6]). The vectors appear to be systematically ordered and gradually curved, suggesting a 
clear hourglass shape that is left-handedly tilted by ~ 60-80° and the direction of the global 
magnetic field that is nearly perpendicular to the elongation of the Serpens cloud cor e from 
NW to SE, ~ 150° (e.g., the C^^O maps of iMcMulhn et al.lboool : loimi fc Testll2002l ). 



Signs of hourglass shapes in the magnet ic field have be en already reported in high- 



mass star forming cores such as NGC 2024 (ILai et al.l 120021). OMC-1 ( ISchleuning I Il998 



Houde et al. 2004: Kusakabe et al. 



2008h. and D R21 Main (|KirbJl2009|). In low mas s cores 



such as NGC 1333 IRAS 4A flGirart et al.ll2006r ) and Barnard 68 flKandori et al.ll2009[ ) hour- 
glass shapes have been more clearly shown. These examples, except OMC-1, of the hourglass- 
shaped magnetic field have been found only in isolated cores or cores with simple structures 
in the star forming regions. However the Serpens cloud core is a molec ular cloud comple x 



consisting of many molecular cloud cores or sub-millimeter cores (e.g., iDavis et al.lll999l ). 
which form a cluster of low-mass YSOs, and the hourglass-shaped magnetic field spreads 
widely over the Serpens cloud core. Thus, this is a clear example that the hourglass-shaped 
magnetic field is associated with a low-mass star forming complex, while OMC-1 is an ex- 
ample of a hight-mass star forming complex. 
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4. Analysis and Discussion 
4.1. Shape of the magnetic field 

We have modeled the shap e of the magn e tic fie ld w ith the polariza t ion ve ctors measured 



at H for point sources, following lGirart et al.l (|2006[ ) and lKandori et al.l (|2009[ ) . The magnetic 
field was fitted with a parabolic function of x = g+gCy"^, with a counterclockwise tilted y-axis 
(the parabolic magnetic field axis of symmetry) by 6'pa and a symmetric center (x, ?/)center, 
where the y is the distance from the horizontal axis {x = 0) and the x is the distance from 
the parabolic magnetic field axis of symmetry. The value of tan'^ {dy / dx) + 90° corresponds 
to the position angle of the polarization (6*). Only the point sources, except YSOs, having 
P/AP > 3 and P < 6.2{H — Ks), were used for the fitting. The error of the polarization 
angle {A6) was used to compute a weight for the datum, 1/(Z\6')^. 

In Figure U\ the best-fit magnetic field is shown as well as the measured polarization 
vectors for 149 sources. The position angle of the parabolic magnetic field axis of symmetry 
is ~ 70°, and the coefiicient C of j/^ is ~ 7.1 x 10~® pixel"^. The root mean square (r.m.s.) 
of the residuals is ~ 22°. 

We executed one-parameter fitting of the magnetic field in local areas, in order to more 
accurately calculate the r.m.s. of the residuals, with the same 6'pa and (x, ?/)center obtained 
in the global fitting above. We selected three corners and one more area of the image 
where the source density is relatively high and/or the magnetic field seems to be rather 
ordered (areas outlined by dashed boxes in Figure [7]). Toward the SE corner of the image 
( X < 400 and y < 400 in Figure O 30 sources), the coefiicient C of y'^ was determined to 
be (7.99 ± 0.76) x 10~^ pixel"^, similar to that of the global fitting, and the r.m.s. of the 
residual was calculated to be 12.9±0.9°, and toward the SW corner (x > 500 and y < 230; 20 
sources), C = (7.52 ±1.00) x 10~^ pixel"^ and r.m.s. = 27.0 ±2.0° were obtained. Removing 
the dispersion due to the measurement uncertainties of the polarization angles 4.2 ±3.0° and 
3.2 ± 2.2°, we obtained the dispersions from the best-fit model, 12.2 ± 1.4°and 26.8 ± 2.0° 
for the SE and SW corners, respectively. Toward the NE corner (x < 300 and y > 800; 18 
sources), C = (3.36 ± 0.92) x 10^^ pixel"^ and r.m.s. = 14.8 ± 1.6° were evaluated, and 
the intrinsic dispersion from our model of 13.7 ± 2.0° was obtained with the measurement 
uncertainty of 5.6 ± 2.4°. This smaller C indicates that the curvature of the magnetic field 
here is rather looser than that expected from the global fitting, i.e., slightly bended to the 
direction parallel to the symmetry axis of the magnetic field. Toward the area next to the 
NE corner (400 < x < 700 and y > 800; 17 sources), C = (6.85 ± 0.55) x 10"^ pixel'^ and 
r.m.s. = 13.3 ± 1.4° were evaluated, and the intrinsic dispersion of 12.6 ± 1.8° was obtained 
with the measurement uncertainty of 4.2 ± 2.9°. 
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4.2. Comparison of the magnetic field with the submilhmeter and miUimeter 

data 



4.2.1. 850 fim continuum 

We compare our H-hand measured polarization vector s and the modeled magnetic field 
with the 850 fim dust continuum map of iDavis et al.l ( 119991 ) in Figure El Note that the green 



lines of this figure do not present lines of magnetic force, just the direction of the magnetic 
field. 

The high intensity ridge of the 850 fim continuum is elongated along the NW-SE di- 
rection, having two sub-clumps (NW and SE sub-clumps), both of which consist of several 
dense cores (e.g., SMM 1-11, S68Nb-d, and PS2 in Figure E]). This distribution of the 
850 fim continuum is very similar to that of the bright parts of the '^'^CQfJ = 1 — 0) and 
C^^O(J =1-0) emission jMcMulhn et al.lboool : loimi fc Testill2002h . ahhough the global 
distributio n of the ^^CO(J = 1 — 0) emission is not always elongated, but rather roundly 
extended (jOlmi &: Testill2002l ). It is evident that the symmetric axis (y'-axis) of the best-fit 
magnetic field with a parabolic function is nearly perpendicular to the elongation direction of 
these continuum and molecular line emissions. The horizontal axis (x'-axis) of the parabolic 
magnetic field is situated nearly along the 850 fim continuum ridge, although there are some 
deviations of the continuum emission from the horizontal axis. The symmetric axis of the 
parabolic magnetic field runs through the northern part of the SE sub-clump, not through 
the middle point of the two sub-clumps, which looks like the center of gravity of the Serpens 
cloud core when we glance at the 850 /xm continuum map. 



Davis et al.l (119991 ) suggested the presence of three extended cavity-like structures to the 



east of SMM 3 (hereafter CLS 1), south-west of SMM 2 (hereafter CLS 2), and north-west 
of SMM 4 (hereafter CLS 3), which consist of three pairs of filaments that protrude the 850 
fim continuum ridge. They mentioned that these cavity structures (CLS 1-3) are probably 
shaped by outfiows rather than by global cloud collapse along, say, magnetic field lines. 

As is in Figure [8], the filaments to the north-east of SMM 3 and east of SMM 2 form 
CLS 1, those to the south-east of SMM2/PS2 and south of SMMll form CLS 2, and those 
to west of SMM3 and east of SMM4 form CLS 3. It appears that the two filaments of CLS 
1 jut almost along the magnetic field from the SE sub-clump and that the symmetry axis 
(j/'-axis) of the magnetic field go through the inside of CLS 1 as well as CLS 3. 
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4-2.2. CO emission 



Here we compare our best-fit magnetic field with 
i^CO J = 1 - 0, and C^^O J = 1 - observations ( 


the ^^CO J = 2 - 


Whi 


te et al. 


1995 


McMullin et al. 


2000; 


Narayanan et al. 


2002; 


Olmi & Testi 


2002). 



1-0, 



12 



CO J = 2 - 1 



As was mentioned above, the bright parts of the ^^CO J = 1 — and C^^O J = 1 — emission 
maps are elongated and confined in the ridge, while the global distribution s of ^^CO J = 2 — 1 



Davis et al. 


1999; 


Olmi & Testi 


2002) 



(e.g. 



White et al. 1995 



Figure [9]present s our best-fit magne tic field superposed on the CO J = 2 — 1 contour map 
and 850 fim image of iDavis et al.l ( 1l999l ). where the CO J = 2 — 1 map is consid ered to show 
the am bient molecular gas of the Serpens cloud core, but not the dense cores. iDavis et al. 
( 119991 ) mentioned that toward the two filaments of CLS 1 and one CLS 2 filament to the 
south-east of SMM2/PS2 the CO J = 2 — 1 emission and 850 /xm continuum distributions 
coincide well. As mentioned above, the two filaments seem to run almost along the magnetic 
field, indicating that the CO J = 2 — 1 filaments are also related with the magnetic field. 
For the CLS 2 filament to the south-east of SMM2/PS2, the same situation as the CLS 1 
filaments may be also seen. Two other CO J = 2 — 1 filaments/extensions to the north-west 
of SMM 9 and west of SMM 1 are also noticeable in Figure O Although considerable parts 
of these two filaments are out of our polarimetry image, the extrapolation of our best-fit 
magnetic field cloud predict that these two filaments run along the magnetic lines. 



C^SQ J =1-0 and^^CO J =1-0 



McMullin et al.l ( l2000l ) showed that a velocity gradient running from a LSR velocity centroid 
of 9 km s~^ at the no rth-west end of the C^ ^O J = 1 — emission to 7.5 km s~^ at the south- 
east end (Figure 2 of iMcMuUin et al.ll2000[ ). i.e., along the elongation direction of C^^O. On 



the other hand, lOlmi Sz Testi f l2002h suggested that the Serpens cloud exhibits a velocity 
gradient roughly from east to west, based on their model fitting of velocity gradients in C^^O 
J = 1 — 0, ^^CO J = 1 — 0, C^'^S J = 1 — 0, adopting their map center, which is the middle 
point of the two sub-clumps, as the reference position for a nalysis. However, a ccording to 
their channel and centroid velocity maps (Figures 7 and 8 of Olmi fc Testi 20021). the brigh t 
parts of C^^O J = 1 - and ^^CO J = 1 - are similar to that of iMcMuUin et al.l fcoooh . 



and a steep velocity gradient from NW to SE almost along the normal line of the symmetry 
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axis (y'-axis) of the magnetic field can be seen at just south of their reference position in 
^^CO, although at the reference position a velocity gradient from West to East is seen. It 
is surprising that the normal line of the steep velocity gradient almost coincides with the 
symmetry axis (y'-axis) of the magnetic field. 

In summary, the direction of velocity gradient is nearly along the elongation of the 
Serpens cloud core and is nearly perpendicular to the symmetry axis of the magnetic field 
with a coincidence of the normal line of the steep velocity gradient and the axis of the 
magnetic field. It could be possible that this normal line of the velocity gradient is an axis 
of the global rotation of the Serpens cloud core if the real center of gravity of the Serpens 
cloud core is located on the symmetry axis of the magnetic field. 

It is interesting to examine the presence of C^^O J = 1 — and ^^CO J = 1 — features 
that coincide with the filaments of the CO J = 2 — 1 e mission and 850 /xm emission . In 
the C^^O 1 - integrated emission maps of I White et al.l ( 1l995l ). iMcMuUin et al.l (120001 ) and 



Olmi fc Testil (120021 ). a feature to the north-east of SMM 3 could coincide with one of the 
CSL 1 filament, but one to the east of S MM 2 is not clear. In the channel map of ^^CO 
J = 1 — (Figure 7 of lOlmi fc Testil |2002| ). a filament feature to the east of SMM 2 is clearly 
visible in the blue-shifted emission at the panel of Vlsr=5-7.3 km s~^. This filament looks 
likely to coincide with the CSL 1 filament to the east of SMM 2, but we can clearly recognize 
that it is located just outside this CSL 1 filament, i.e., between this CSL 1 filament and 
the CSL 2 filament to the south-east of SMM 2/PS2. At the same panel, a feature to the 
north and north-east of SMM 3 or near SMM 8 is also visible. This feature appears to be 
just outside the CLS 1 filament to the north-east of SMM 3. At the panel of Vlsr=S.4:- 
12 A km s~^, a red-shifted feature that protrudes from the SE sub-clump is visible, but it is 
located toward the inside region of CLS 1. The presence of this red-shifted feature and the 
blue-shifted features are probably consistent with red-shifted velocity region that jut from 
the SE sub-clump and with blue-shifted regions toward both sides of this r e d-shif ted region, 
respectively, in the ^^CO J = 1 — centroid velocity map of lOlmi &: Testil (120021 ). 



CO outflows 



Davis et al.l ( 1l999l ) presented the integ rated intensity con tours of CO J = 2 — 1 blue- and 



red-shifted outflows (Figures 4 and 8 of iDavis et al.lll999l ). These figures imply that the 850 
ixm. filaments that coincide the CO J = 2 — 1 filaments are shaped by outflows. On the basis 
of a fact that these filaments run along the magnetic field, the outflows that protruded from 
the ridge to its ambient are most likely to be guided by the magnetic field or to drag the 
magnetic field. The outflows may be guided by the magnetic field since the magnetic field 
seems to be strong enough to be ordered at least over our polarimetric imaging area. 
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J 



The CLS 1 filaments are associated with red-shifted outflows, but no red-shifted CO 
2 — 1 emission is visible at the root of CLS 1. However, CO J = 1 — obervations 



(INarayanan et al.ll2002l ) showed U-shaped, red-shifted high velocity flow at the root of CLS 
1. Th is CO J = 1 — feature and our best-fit magnetic field support the idea of iDavis et al. 
(119991 ) that the CLS 1 filaments of the CO J = 2 — 1 and 850 /im emission illustrate the 
action of a wide-angled wind powered by a source within the SMM 2/3/4 cluster, which has 
swept up gas and dust into a warm, compressed shell, although there is a possibility that 
the wind is powered by multiple sources within the cluster. 



4.3. Magnetic field strength 



We try to make an evaluation of the magnetic field strength toward four areas where we 
calculated the angular dispersions (residuals) for our best-fit mag netic field, using the Chan- 



drasekhar & Fermi (CF) method (jChandrasekhar fc Fermilll953l ). On the basis of the con- 



clusions of recent MHD studies that the introduction of a correction factor is needed for eval- 



(Ostriker et al. 


2001 




Padoan et al. 


(2004 


) mentioned that a correction 



factor of ~0.5 is appropriate in most cases when the magnetic field is not too weak. Since 
the magnetic field seems to be ordered over the Serpens cloud core, the magnetic field is 
expected to be strong. Therefore we first adopt a correction factor of 0.5 to evaluate the 
magnetic field strength. We need the mass density and velocity dispersion of the matter 
coupled to the magnetic field to eval uate the magrietic fi eld strength. Here, we use those 



estimated from the C^^O observation (jOlmi fc Testill2002l ). 



Toward the four areas, the Ho column densities from C O could be estimated to be 



6 X 10^^ cm ^ from Figure 11 of lOlmi fc Testil (|2002f). Adopting the approximate C^^O 



extent of ~12'(~0.9 pc at d~260 pc; Figure 2 of lOlmi k Testil l2002f ) as the depth of thes e 
area, we obtain the H2 densities of ~ 2.1 x 10^ cm~^. From Figure 10 of lOlmi &: Testil (120021 ) . 
the C^^O velocity widths could be estimated to be ~1.6-1.8 km s^^ toward the SE and NE 
corners, and ~1. 8-2.0 km s~^ toward the area next to the NE corner. Toward the SW corner 
with a complex distribution of velocity width, the velocity width may be ~l-2 km s~^. 
Using a mean molecular mass, /i, of 2.3 and these values to derive the velocity dispersions, 
we roughly evaluated the magnetic field strength of the plane-of-the-sky of B\\ ~ 160-180 
jiG toward the SE corner, ~ 150-160 /xG toward the NE corner, and ~ 180-200 /zG toward 
the area next to the NE corner. Although ~50-90 ixG can be evaluated toward the SW 
corner, this value might be more uncertain than those toward the other areas due to the 
larger uncertainty of the velocity width. The magnetic field strength evaluated here is 
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highe r than those measured around dark cl oud complexes and prestellar cores, a few 10 /iG 



Alves et al. 2008: Kandori et a 



20091 respectively), but smaller than thos e around HII 



(e.g., 

regions, a few mG (e.g.. lHoude II2004I ) and of a protostellar envelope, a few mG (iGirart et al. 

mm- 



Recently, iHoude et al.l (|2009[ ) showed how the signal integration through the thickness of 
the cloud and the area of the telescope beam affects on the measured angular dispersion and 
apply their results to OMC-1. Based on their estimated number (N=21) of the independent 
turbulent cells contained within the column probed by the telescope beam, they found that 
a correction factor of 1/\/N ~ 0.2 is applicable to OMC-1. In our case, although the area 
of the telescope beam is negligibly small due to the point sources, the thickness of the cloud 
should be taken into account and the correction factor should be somewhat smaller than 
~ 0.5. If we assume that the effect of the cloud thickness is similar to that of OMC-1, 
we obtain ~ 11, suggesting a factor of ~ 0.3. Adopting this factor of ~ 0.3, the above 
estimated values are reduced by a factor of ~ 0.6 and B\\ ~ 100 may be appropriate for the 
ambient region of the Serpens cloud core, except the SW corner. 

Here, we roughly derive the mass to magnetic flux ratio Mdoud/^ using our estimated 
value of B ~100 /iG, and compare it wit h the critical value for a rn agnetic stability of 
the cloud, (Mcioud/^)criticai = (47r2G)-i/2 flNakano fc Nakamura I llQTsI ). With a formula 
A^cioud/^ = {nR'^fimuN)/{7iR'^B) = finiuN/B and the H2 column density A^ ~ 6 x 10^^ 
cm~^ where we estimated B, we derive Mdoud/^ ~ 3.8 x (Mdoud/^) critical, where i? is a 
radius of the cloud and mn is the mass of a hydrogen atom. Although this derived value is 
slightly larger than the critical value, Mdoud/^ could be much larger in the inner region of the 
cloud core because the column density of the inner region is much higher than those where 
we estimated B, but the magnetic field may be slightly larger than that we estimated in the 
ambient region, judged from the slowly curved shape of the magnetic field. We note that the 
adopted strength of the magnetic field is that estimated for the projection of the magnetic 
field in the plane of the sky, suggesting a slightly smaller Mdoud/^ than the estimated one. 
These imply that the ambient region is marginally supercritical, while the inner region is 
supercritical. This situation is considered to be quite consistent with the hourglass shape of 
the magnetic field and with the cluster formation within the sub-clumps. 

It is interesting to examine whether the magnetic field can maintain the outflow colli- 
mation along the magnetic field in the ambient region of the sub-clumps, i.e., whether the 
magnetic field can guide the outflows. The magnetic pressure, Pg = B'^/8n, is calculated to 
be ~ 4 xlO~^° dyn, adopting B ~ 100/iG. Assuming the average density and velocity width 
due to turbulence for the outflow to be 3 x 10^ cm~'^, which woul d be consistent with the 
optically thin condition of the high velocity gas (IWhite et al.lll995l ). and 3 km s~^, which is 
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larger than the C^^O velocity width by a factor of ~1. 5-2.0, we obtain the turbulent pressure, 
-Pturb = pcTturb) ~ 2xl0~^° dyu. Taking into account the fact that the adopted strength 
of the magnetic field is that estimated for the projection and that is proportional to B^, 
these estimates imply that the magnetic field can guide the outflows in the ambient region 
of the Serpen cloud core. 



4.4. Comparison with outflow-driven turbulence model for cluster formation 



From our analysis, the magnetic field seems to be important in considering the cloud 
stability that is related to star formation or cluster formation and the feed back from the 
star formation activity, such as outflows. 

The hourglass-shaped magnetic field suggests that the Serpens cloud core first contracted 
along the straight magnetic field to be a filament or elongated cloud, which is perpendicular 
to the magnetic field, and that then the central part contracted cross the magnetic field due 
to the high density in the central region of the cloud core. This situation is very similar to 
the contraction of the low- mass core t hat is penetrated by the uniform magnetic field (e,g., 
Girart et al.ll2006l : iKandori et al.ll2009l ). In addition, there might exist the cloud rotation, of 
which axis agrees with that of the hourglass-shaped magnetic field. It was report ed that many 



small-scale outflows spread to or penetrate the NW and SW sub-clupms(e.g., iHerbst et al. 



19971 : iHodapp Ill999l : iDavis et al.lll999l : IZiener fc Eisloffellll999l ) , and the ar abient, larger-scale 



outflows (filamerits) se em to run along the magnetic field as shown above (IDavis et al.lll999 
Narayanan et al.ll2002l ). Moreover, it is possible that the blue-shifted ^'^CO (1 — 0) features 
just outside CLS 1, which correspond to the red-shifted CO (2 — 1) outflows, are inflows 
from the ambient to the central part of the SE sub-clump. Considering these altogether, 
we may have to take into account the magnetic field, outflows, inflows, cloud rotation, and 
contraction as well as the turbulence of the molecular gas in the cluster formation process 
of the Serpens cloud core (see Figure [TOj) . 

The structures mentioned above seern . to be in good agreeme r it wit h the outflow-driven 
turbulence modelof iLi fc Nakamural (120061 ) and iNakamura fc Li I (120071 ) who performed 3D 
MHD simulation of cluster formation takinginto accout the effect of protostellar outflows. 
They demonstrated that protostellar outflows can generate supersonic turbulence in pc-scale 
cluster forming clumps like the Serpens cloud core. One of the important characteristics 
of outflow-driven turbulence is that gravitational infall motions almost balance the outward 
motions dri ven by outflows, creatin g very comphcated density and velocity structure (see e.g.. 
Figure 4 of iNakamura &: Li II2007I ). The resulting quasi-equilibrium state can be maintained 
through active star formation in the central dense region. In the presence of relatively 
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strong magnetic field, both outflow and inflow motions in the less dense envelope tend to 
be guided by large scale ordered magnetic field lines. As a result, filamentary strucutures 
that are roughly converging toward the central dense region appear in the envelope, whereas 
the density structure tends to be more complicated in the central dense region where self- 
gravity and turbulence may dominate over the magnetic field. Infall motions detected by 
^^CO (1 — 0) in the Serpens core may correspond to such filamentary structures created by 
gravitational infall. 

To clarify how the outflows and magnetic fi eld affect the dynarn ical state of the cloud, 
we assess the force balance in the cloud, following iMaury et al.l (120091 ). To prevent the global 
gravitational contraction, the following pressure gradient is needed to achieve the hydrostatic 
equilibrium: 



grav 



dr 



^ M{r)p{r) 



a 



-2 



) 



(1) 



where M(r) is the mass contained within the radius r and we assume that the cloud is 
spherical. The effect of magnetic field is taken into account by the factor (1 — and a is 
the mass-to-magnetic flux ratio normalized to the critical value and is approximated as 



a 



27rGV2M/7rr2 
B 



(2) 



[e.g. 



Nakanol 119981 ). 



Assuming the density profile of p oc r ^, the pressure needed to support the cloud 
against the gravity is estimated to be 

G'M(i?)2 



P ~ 

^ grav 



(1 - «-) 



The force needed to balance the gravitational force is thus evaluated to be 

GM{R)^ 



grav — 



(1-a-) . 



(3) 



(4) 



Adopting M{R) = 2IOM0, R = 0.46 pc flOlmi fc Testill2002f ). B = lOOpG and a = 3.8, Fg^^v 
can be estimated to be ~ 4.3 x 1O~^M0 km s~^ yr~^. The moderately strong magnetic field 
of a = 3.8 can reduce the gra vitational forc e by ~ 7%. We note that we rescaled the cloud 
mass and radius derived from lOlmi fc Testil (120021 ) by assumingthe distance to the cloud of 
260 pc. Hereafter, we also use other values rescaled for this distance. 



On the basis of the CO (J = 2 — 1) observations, iDavis et al.l (119991 ) detected many 
powerful CO outflows in this cloud, and derived the physical properties of the outflows. From 
their analysis, we can evaluate the total force exerted by the outflows in this region as 



outflow 



Poutflow 8.7-17.5M0km s" 



2.5 X lOV 



(3.4-7.0) X IO-^Mq km s"^ yr^^ 



(5) 
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where Poutfiow is the total outflow momentum, and tdyn is the representative djTiamical time 
of the outflows. The force due to the outflows, -Foutflow, is comparable to or somewhat 
larger than the force needed to stop the global gravitational collapse, -Fgrav, suggesting that 
the outflows play a crucial role i n the cloud dynamics. This result, however, apparently 
contradicts that of lOlmi fc Testil (120021 ) who suggested that the cloud may be undergoing a 
global contraction, although the further justification is needed to confirm their interpretation. 
This apparent inconsistency may come from our assumption of the spherical cloud. Since the 
relatively strong magnetic field associated with the cloud can guide the large scale outflow 
motions along the global magnetic field as discussed in the previous subsection, the force 
exerted by the outflows is expected to be weak along the cross- field direction. As a result, 
the cloud may be able to contract along the cross-field direction. For the Serpens core, both 
the magnetic field and the outflows are likely to control the cloud dynamics. 

The outflows are also expected to be the major source for generating supersonic tur- 
bulence in t he Se rpens core. From the physical quantities of the outflows measured by 
Davis et al.l ( 1l999l ). we can evaluate the total energy injection rate due to the outflows in 
this region as 



dE. 



outflow 



outflow 



^12.7-48.3) J 



(0.5-2)Lc 



dt tdyn 2.5 X 10^ yr ' ' ® ' 

where -^outflow is the total outf l ow en ergy. The energy dissipation rate of supersonic turbu- 
lence is obtained by iMac Lowl (Il999[ ) as 



dEturh ,.l/2MAV^ 



dt 



f 



(7) 



where /(= 0.34) is the non-dimensional constant determined from the numerical simulations, 
and M is the cloud mass, and AV^ is the ID FWHM velocity width. The driv ing scale of the 
turbulence is estima ted to be ~ 0.4 pc for the outflow-driven turbulence (lMatzneiil2007l : 



Nakamura &: Li 1120071 ). The energy dissipation rate of the turbulence can be estimated to 
be ~ O.I2L0, where the FWHM velocity width of about 2 km s"^ is adopted (jOlmi &: Testi 



20021 ). This energy dissipation rate is somewhat smaller than the outflow energy input rate. 
In the Serpens cloud core, the relatively strong magnetic field tends to guide the outflows 
and therefore the significant amount of the outflow energy might escape away from the 
cloud along the magnetic field, as inferred from the magnetic field and outflow structures 
discussed above. In any case, the outflows seem to have sufficient energy to power supersonic 
turbulence in this region and the magnetic field seems to play an important role in the escape 
of the outfiow energy from the cloud. These characteristics appear to be in agreement with 
the the outfiow-driven turbulence model for cluster formation, and imply the importance of 
the magnetic field for the continuous star formation in the central region of the Serpens cloud 
core under the condition where the outfiow energy injection rate is high. The Serpens cloud 
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core is expected to be one of the good examples of the outflow-driven turbulence model for 
cluster formation. 



4.5. Summary 

We have conducted deep and wide (~7f7 x T.l) JHKg imaging polarimetry of the 
Serpens cloud core. The main findings are as follows: 

1. The central part of the infrared reflection nebula is illuminated mainly by two sources; 
the north by SVS 2 (SRN) and the south by SVS 20 with two centrosymmetric patterns. 
The characteristics of the nebula are consistent with those reported in the previous infrared 
polarimetric works. Detailed inspection enabled us to find 24 YSOs associated with IR 
nebulae, in addition to SVS 2 and SVS 20. 

2. Polarization of NIR point sources was measured and those sources, except YSOs, 
have an upper limit of polarization degree similar to that of the nearby star forming regions. 
It is consistent with the dichroic origin, i.e., the polarization vectors of the near-IR point 
sources could indicate the direction of the averaged local magnetic field. 

3. The polarization vectors suggest a clear hourglass shape. We have made a model 
fitting of this shape with a parabolic function and found that the symmetry axis (^pa ~70°) 
of the hourglass magnetic field is nearly perpendicular to the elongation (~ 150°) of the 

bright parts of C^^O ( J = 1 — 0) or submillimeter continuum emissions, i.e., the alignment 
direction of NW and SE sub-clumps. The submillimeter continuum filaments and CO outflow 
lobes, which protrude from these sub-clumps, seems to run along the best-fit magnetic field 
in the ambient region and some ^^CO velocity features also seem to be along the magnetic 
field. 

4. The evaluation of the magnetic field strength has been done with the CF method 
toward the ambient area of the Serpens cloud core, taking into account the recent study on 
the signal integration effect for the dispersion component of the magnetic fleld. The mass 
to magnetic flux ratio was estimated with the evaluated magnetic fleld strength of ~ lOO/xG 
and the parameters of the previous C^^O ( J = 1 — 0) observations, and found to be slightly 
larger than the critical value of magnetic instability in the the ambient area. This suggests 
a possibility that the central region is magnetically unstable, which is consistent with the 
fact that star formation is actively taking place in the central region. We estimated the 
magnetic pressure and the turbulent pressure of the outflow using the evaluated magnetic 
fleld strength and possible turbulent parameters, and found that the magnetic pressure could 
be high enough to guide the outflows in the ambient region. 
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5. The bright part of C^^O (J = 1 — 0), submilhmeter continuum cores as well as 
many class O/I objects are located just toward the constriction region of the hourglass- 
shaped magnetic field. These suggest that the Serpens cloud core first contracted along the 
magnetic field to be an elongated cloud and that then the central part contracted cross the 
magnetic field due to the high density in the central region of the cloud core. 

6. Comparisons of the best-fit magnetic field with the previous observations of molecular 
gas and large-scale outflows suggest a possibility that the cloud dynamics is controlled by 
the magnetic field, protostellar outflows and gravitational inflows. In addition, the outflow 
energy injection rate appears to be the same as or larger than the dissipation rate of the 
turbulent energy in this cloud, indicating that the outflows are the main source of turbulence 
and that the magnetic field plays an important role both in allowing the outflow energy 
to escape from the central region of the cloud core and enabling the gravitational inflows 
from the ambient region to the central region. These characteristics appear to be in good 
agreement with the outflow-driven turbulence model for cluster formation and imply the 
importance of the magnetic fleld to continuous star formation in the center region. 

We thank Chris Davis for kindly providing us the 850 fim continuum data. This work 
was supported by Grant-in-Aid for Scientiflc Research (20403003, 19540242, and 19204018) 
from the Ministry of Education, Culture, Sports, Science and Technology. 

A. Identification of YSOs with NIR reflection nebulae toward the central 

region of the Serpens cloud core 

Except isolated YSOs that are on the periphery of the central region, it is not easy 
to examine whether YSOs have reflection nebulae locally illuminated by themselves only 
with the binned map of the K s polarization vectors (Figure [3]), due to the contamination 
of light from the strong sources or nearby sources. We constructed the highest resolution 
maps without binning for the sources suffering from the contamination of SVS 2, SVS 20 and 
the members of the SVS 4 cluster, and tried to identify which sources are associated with 
reflection nebulae. In Figure [TT| we show the polarization vector maps only for YSOs that we 
have identifled as those having reflection nebulae. For EC 94, EC 98, and EC 121, although 
the vector map quality /resolution is not always good enough for the robust identiflcation, 
we concluded, taking account the weak emission around these sources, that they probably 
have reflection nebulae. 
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Fig. 1. — DSS II R image of the Serpens cloud core. The area of the polarimetric imaging 
is indicated by a rectangle with a size of 8.4' x 8.5' and a center position of {a, 5)j2ooo 
=(18:29:57.6, +01:14:34). 



Fig. 2. — JHKs polarization vector maps of the Serpens cloud core, superposed on the total 
and polarized intensity images in the logarithmic scale. The vectors were made by 12x12 
pixel binning. The field of view is ~ 7 '.7 x 7'.7. North is at the top and east is to the left. 
The center position of each image is (q;,5)j2ooo =(18:29:57.39, +01:14:36.2). The dead pixel 
regions of the J-band images are at the upper right and near the middle of the right edge. 
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Fig. 3. — Ks polarization vector map of the central region of the Serpens cloud core, su- 
perposed on the total intensity images in the logarithmic scale. The vectors were made by 
3x3 pixel binning. The area of the image is 220" x 220". The reference position of offset is 
(a, ^)j2ooo =(18:29:57.39, +01:14:36.2). North is at the top and east is to the left. 40% vector 
is shown near the bottom right corner. The young stellar objects with reflection nebulae are 
marked. The SVS 4 cluster is shown by the circle. 
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Fig. 4. — Ks polarization vector map of the NW sub-core region of the Serpens cloud core, 
superposed on the total intensity images in the logarithmic scale. The vectors were shown 
every 2 pixels. The area of the image is 130" x 130". The reference position of offset is 
(a,'^)j2ooo =(18:29:57.39, +01:14:36.2). North is at the top and east is to the left. 40% 
vector is shown near the bottom right corner. The young stellar objects with reflection 
nebulae are marked, but EC28 and EC41 are not identified as young stellar objects with 
reflection nebulae. 
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Fig. 5. — Top: polarization degrees at H versus H — Ks color for sources having the H-hand 
polarization errors of < 0.3%. Bottom: polarization P.A.s at H versus H — K color for 
sources having P/ AP > 3 and P < 6.2{H — Ks). 



Fig. 6. — //-band polarization vectors (white bars) for point sources having Pj AP > 3 
and P < 6. 2 (if — Ks), superposed on a JHKs composite color image in the logarithmic 
scale. 10% vector is shown near the upper right corner. The area of the image (~ 8'A x 8^5) 
made by averaging the dithered frames is slightly larger than that of the camera array 
(~ 7'.7 X 7'.7). North is at the top and east is to the left. YSOs identified by the Spitzer space 
telescope are marked by circles (Class O/I; red, Flat-spectrum; magenta, C lass II; yellow, 



Transition disk; blue. Class III; white) with ID numbers (Tables 3 and 4 of IWinston et al. 



20071 ). Submillimeter continuum peaks are show by cros s es wit h names and their coordinates 



come from the Sp itzer photometry of IWinston et al.l (120071 1 . except SMM2 and SMMll 
flPavis et al1ll999[ ). 
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Fig. 7. — if-band measured polarization vectors (solid bars) for point sources having 
P/AP > 3 and P < 6.2{H — Ks) and their best-fit polarization (magnetic field) vec- 
tors (dashed bars) with a parabolic function {x = g + gCy"^), where the x is the dis- 
tance from the parabolic magnetic field axis of symmetry having the coefficient of y'^ 
oi C = (7.00 ± 0.39) X 10~^ pixel"^ and the symmetry center of the magnetic field is 
{x,y)center = (646.0 ± 11 .3, 384.0 ± 9. 1) , which corresponds to (a,(5)j2ooo =(18:29:54.9, 
+01:13:13), on the axis of symmetry that is tilted at 9pA = 70.01° ± 0.77°. 10% vector 
is shown near the upper right corner. The area of the figure is ~ 8'A x 8^5, nearly same as 
that of Fig. O YSOs identified by the Spitzer space telescope are marked by circles, and 
these YSOs are not used for the fitting. In four areas enclosed by dashed lines, additional 
fitting with one-parameter (C) was done (see text). 



Fig. 8. — //-band polarization vectors (white lines) for point sources having Pj AP > 3 
and P < 6.2{H — Ks), and their best-fit magnetic field (green thick curved/ straight lines 
and thin curved lines), super posed on SCUBA 850 /xm continuum image (kindly provided 
by C.J. Davis, see Figure 1 of iDavis et al.l 119991 ). Note that these lines do not present lines 
of magnetic force, just the direction of the magnetic field. The horizontal axis (green thin 
straight line), which is perpendicular to the parabolic magnetic field axis (green thick straight 
line), is also shown. The area of the image is ~ 8^4 x 8^5, same as that of Fig. |6l North is 
at the top and east is to the left. The YSOs identified by the Spitzer space telescope and 
submillimeter continuum peaks are indicated the same as in Figure Ei 
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Fig. 9. — The best-fit magnetic field (the same as Figure [8]), superposed on the CO J = 2 — 1 
conto ur plot of Vlsr=2-16 km s~^ and 850 fim continuum image (Figure 2 of iDavis et al. 
19991 ). 




Fig. 10. — The schematic drawing of the Serpens cloud core. The pinched or hourglass- 
shaped magnetic field pattern suggests that the magnetic field lines are dragged along with 
the contracting gas toward the center of the cluster. 



Fig. 11. — Ks polarization vector maps of the nebulosities associated with YSOs in the SVS 
4 cluster and the central region of the Serpens cloud core, superposed on the total intensity 
images in the logarithmic scale. The vectors were shown every 2 pixels. The reference 
position of offset is («,5)j2ooo =(18:29:57.39, +01:14:36.2). North is at the top and east is 
to the left. 40% vector is shown at each panel. YSOs with reflection nebulae are marked. 
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Table 1. YSOs with NIR nebulae toward the central region. 



YSO Name'' 




Source 


SSTc2dJ 


ISO 


Other Name'' 


(Spitzer ID'' ) 


(mag)'' 


Class^ 


ID° 


IDd 




EC79 (80)^ 


11.5 


2 


18295655+0112595 


304 


GCNM84/STGM12 


SVS2 (9) 


9.3 


0/1 


18295687+0114465 


307 


EC82/GCNM87/CK3/STGM22 


EC84 (85)° 


11.1 


2 


18295696+0112477 


309 


GCNM90/STGM11 


EC89 (12)° 


12.0 


0/1 


18295766+0113046 


(312) 


GCNM97/STGM13 


SVS20S/N (35) 


7.1 


FS 


18295772+0114057 


314 


EC90/GCNM98/CK1/STGM18 


EC91 (70)" 


13.0 


2 


18295780+0112279 


320 


GCNMIOI 


EC92 (2)« 


10.5 


0/1 


18295783+0112514 


(317) 


GCNM104 


EC93 (83) 


10.8 


2 


18295780+0115318 


319 


GCNM100/STGM25/CK12 


EC94 (37) " 


11.7 


FS 


18295784+0112378 


318 


GCNM102 


EC95 (105)'= 


10.0 


2 


18295789+0112462 


(317) 


GCNM103 


EC97 (27) 


9.9 


FS 


18295819+0115218 


321 


GCNM106/CK4/STGM24 


EC98 (165)° 


12.3 


TD 


18295844+0112501 


322 


GCNMllO 


EC103 (4) 


11.8 


0/1 


18295877+0114262 


326 


GCNM112/STGM20/K2.5 


EC105 (59) 


9.5 


2 


18295923+0114077 


328 


GGNM119/CK8/STGM19/K2.6 


EC117 (216) 


10.1 


3 


18300065+0113402 


338 


GCNM135/CK6 


EC118 (158) 


9.0 


TD 


18300061+0115204 


337 


GCNM136/CK2 


EC121 (30) 


13.2 


FS 


18300109+0113244 


341 


GCNM142 


EC125 (1) 


13.4 


0/1 


18300208+0113589 


345 


GCNM154/CK7/STGM16 


EC129 (10) 


9.9 


0/1 


18300273+0112282 


347 


GCNM160/STGM10 


''Prom the names 


referred 


in Table 1 


of lEiroa et al.l ||200S|~1 







''From lWinston et all fcOOTl) 
°From rHarvev et al ] l|2007h 
'^Prom lKaas et all l|2004l ) 



°SVS-4 cluster 
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Table 2. YSOs with NIR nebulae toward the NW region. 



YSO Name'' 


Ks 


Source 


SSTc2dJ 


ISO 


Other Name'' 


(Spitzer ID'' ) 


(mag)'' 


Class'' 


ID" 


IDd 




DEOS/S68Nd<= (11) 


14.9 


0/1 


18294913+0116198 


250 


knot c' /K4.5/WMW11 


S68Nc° 




(0/1) 






subknot a3'' /knot a^ 


SMMl-FIRSl 






18294963+0115219 


258a 


GCNM23 /WMW114/VLA7 


EC53/SMM5 (24) 


11. S'^ 


0/1 


18295114+0116406 


265 


STGM27/WMW24 


EC67 (81) 


9.6 


2 


18295359+0117018 


283 


GCNM60/STGM29/WMW81 


EC38/S68Nb° (7) 


12.7 


0/1 


18294957+0117060 


254 


WMW7 


SMMIO-IR (21) 


17.7^^ 


0/1 


18295219+0115478 


270 


WMW21 



''From the names referred in Table 1 of lEiroa et al. 1 1 I2OO8I) 

faProm lWinston et al.l ||2007|) 

'^From lHarvev et al.l ll2007h 

''Prom lKaas et al.l ll2004h 

"From lWilliams fc Mver3 ll2000h 

^From lPavis et al. I lll999l) 

gFrom lHodappI lll999l) 



